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Abstract 

Using the characterization of last multipliers as solutions of the 
Liouville's transport equation, new results are given in this approach 
of ODE by providing several new characterizations, e.g. in terms of 
Witten and Marsden differentials or adjoint vector field. Applications 
to Hamiltonian vector fields on Poisson manifolds and vector fields on 
Riemannian manifolds are presented. In Poisson case, the unimodular 
bracket considerably simplifies computations while, in the Riemannian 
framework, a Helmholtz type decomposition yields remarkable exam- 
ples: one is the quadratic porous medium equation, the second (the 
autonomous version of the previous) produces harmonic square func- 
tions, while the third refers to the gradient of the distance function 
with respect to a two dimensional rotationally symmetric metric. A fi- 
nal example relates the solutions of Helmholtz (particularly Laplace) 
equation to provide a last multiplier for a gradient vector field. A 
connection of our subject with gas dynamics in Riemannian setting is 
pointed at the end. 
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Introduction 



In January 1838, Joseph Liouville(1809-1882) published a note ([10]) on the 
time-dependence of the Jacobian of the "transformation" exerted by the so- 
lution of an ODE on its initial condition. In modern language, if A = A(x) is 
the vector field corresponding to the given ODE and m = m(t, x) is a smooth 
function (depending also of the time t), then the main equation of the cited 
paper is: 

^ + m ■ divA = {LE) 

called, by then, the Liouville equation. Some authors use the name general- 
ized Liouville equation ([4]) but we prefer to name it the Liouville equation 
of transport (or of continuity). This equation is a main tool in statistical 
mechanics where a solution is called a probability density function ([22]). 

The notion of last multiplier, introduced by Carl Gustav Jacob Jacobi 
(1804-1851) around 1844, was treated in detailed in "Vorlesugen fiber Dy- 
namik", edited by R. F. A. Clebsch in Berlin in 1866. So that, sometimes it 
has been used under the name of " Jacobi multiplier" . Since then, this tool 
for understanding ODE was intensively studied by mathematicians in the 
usual Euclidean space M n , cf. the bibliography of [3], [13]-[16]. For all those 
interested in historical aspects an excellent survey can be found in [2]. 

The aim of the present paper is to show that last multipliers are ex- 
actly the autonomous, i.e. time-independent, solutions of LE and to discuss 
some results of this useful theory extended to differentiable manifolds. Our 
study has been inspired by the results presented in [17] using the calculus on 
manifolds especially the Lie derivative. Since the Poisson and Riemannian 
geometries are the most frequently used frameworks, a Poisson bracket and a 
Riemannian metric are added and cases yielding last multipliers are charac- 
terized in terms of unimodular Poisson brackets and respectively, harmonic 
functions. 

The paper is structured as follows. The first section reviews the defini- 
tion of last multipliers and some previous results. New characterizations in 
terms of de Rham cohomology and other types of differentials than the usual 
exterior derivative, namely Witten and Marsden, are given. There follows 
that the last multipliers are exactly the first integrals of the adjoint vector 
field. For a fixed smooth function m, the set of vector fields admitting m as 
last multiplier is a Lie subalgebra of the Lie algebra of vector fields. 



2 



In the next section the Poisson framework is discussed, on showing that 
some simplifications are possible for the unimodular case. For example, in 
unimodular setting, the set of functions / which are last multiplier for exactly 
the Hamiltonian vector field generated by / is a Poisson subalgebra. Local 
expressions for the main results of this section are provided in terms of the 
bivector n defining the Poisson bracket. Two examples are given: one with 
respect to dimension two when as a last multiplier is obtained exactly the 
function defining the bivector, and the second related to Lie-Poisson brackets 
when the condition to be last multiplier is expressed in terms of structural 
constants of given Lie algebra. Again, the two dimensional case is stud- 
ied in detail and for non-vanishing structure constants, an afline function is 
obtained as a last multiplier. 

The last section, is devoted to the Riemannian manifolds and again, some 
characterizations are given in terms of vanishing of some associated differen- 
tial operators, e.g. the co differential. Assuming a Helmholtz type decomposi- 
tion, several examples are given: the first is related to a parabolic equation of 
porous medium type and the second yields harmonic square functions. As to 
the first example, one should notice the well-known relationship between the 
heat equation (in our case, a slight generalization) and the general method 
of multipliers; see the examples from [19, p. 364]. The third example is de- 
voted to the distance function on a two dimensional rotationally symmetric 
Riemannian manifold, while the final example is concerned with solutions of 
Helmholtz (particularly Laplace) equation, for obtaining a last multiplier for 
a gradient vector field. 

The last section deals with an application of our framework to gas dy- 
namics on Riemannian manifolds but in order to not enlarge our paper too 
much we invite the reader to see details of physical nature in Sibner's paper 
[20]. More precisely, a main result from [20] is rephrased in terms of last 
multipliers. 

Acknowledgments The author expresses his thanks to loan Bucataru, 
Marian-loan Munteanu, Zbiegnew Oziewicz and Izu Vaisman for several use- 
ful remarks. Also, David Bao and an anonymous referee suggest some im- 
portant improvements. 



3 



1 General facts on last multipliers 



Let M be a real, smooth, n-dimensional manifold, C°° (M) the algebra of 
smooth real functions on M, X (M) the Lie algebra of vector fields and 
A fc (M) the C°° (M)-module of ^-differential forms, < k < n. Assume that 
M is orientable with the fixed volume form VeA" (M). 
Let: 

x (t) = A 1 (x l (t),...,x n (t)) , 1 < i < n 

an ODE system on M defined by the vector field A e X (M) , A = (A') 1 < i < n 
and let us consider the (n - l)-form Q = i A V e A™ -1 (M). 

Definition 1.1([6, p. 107], [17, p. 428]) The function m G C°° (M) is 
called a last multiplier of the ODE system generated by A, (/as£ multiplier 
of A, for short) if mVl is closed: 

d(mfi) := (dm) AQ + mdQ = 0. (1.1) 

For example, in dimension 2, the notions of last multiplier and integrating 
factor are identical and Sophus Lie suggests a method to associate a last 
multiplier to every symmetry vector field of A (Theorem 1.1 in [8, p. 752]). 
Lie method is extended to any dimension in [17]. 

If the (n — l) th de Rham cohomology space of M is zero, H n ~ l (M) = 0, 
there follows that m is a last multiplier iff there exists a E A n ~ 2 (M) so that 
mfl = da. Other characterizations can be obtained in terms of Witten's 
differential [25] and Marsden's differential [11, p. 220]. If / e C°° (M) and 
t > 0, Witten deformation of the usual differential d tf : A* (M) -> A* +1 (M) 
is defined by: 

dt f = e~ tf de tf 

which means [25]: 

d t f (u) = tdf Aw + duj. 
Hence, m is a last multiplier if and only if: 

d m fl = (1 — m) dfl 

i.e. fl belongs to the kernel of the differential operator d m + (m — 1) d : 
A^ 1 (M) -> A n (M). Marsden differential is d f : A* (M) -> A* +1 (M) de- 
fined by: 

df {u) = -d{fu) 
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and m is a last multiplier iff Vt is cf"-closed. 

The following characterization of last multipliers will be useful: 
Lemma 1.2([17, p. 428]) (i) m G C°° (M) is a last multiplier for A if 

and only if. 

A{m) + m-divA = (1.2) 

where divA is the divergence of A with respect to volume form V. 
(ii) Let ^ h G C°° (M) such that: 

L A h := A (h) = (divA) ■ h (1.3) 

Then m = h^ 1 is a last multiplier for A. 

Remarks 1.3 (i) Equation (1.2) is exactly the time-independent version 
of LE from the Introduction. So that, the promised relationship between LE 
and last multipliers is obtained. An important feature of equation (1.2) is 
that it does not always admit solutions cf. [7, p. 269]. 

(ii) In the terminology of [2, p. 89], a function h satisfying (1.3) is called 
an inverse multiplier. 

(iii) A first result given by (1.2) is the characterization of last multipliers 
for divergence-free vector fields: m G C°° (M) is a last multiplier for the 
divergenceless vector field A iff m is a first integral of A. The importance of 
this result is shown by the fact that three remarkable classes of divergence-free 
vector fields are provided by: Killing vector fields in Riemannian geometry, 
Hamiltonian vector fields in symplectic geometry and Reeb vector fields in 
contact geometry. Also, there are many equations of mathematical physics 
corresponding to the vector fields without divergence. 

(iv) For the general case, namely A is not divergenceless, there is a strong 
connection between first integrals and last multipliers as well. Namely, from 
properties of Lie derivative, the ratio of two last multipliers is a first integral 
and conversely, the product between a first integral and a last multiplier is 
a last multiplier. So, denoting Flnt(A) the set of first integrals of A, since 
Flnt(A) is a subalgebra in C°°(M) it results that the set of last multipliers 
for A is a F Int(A)-module. 

(v) Recalling formula: 

div(fX)=X(f) + fdivX (1.4) 

there follows that m is a last multiplier for A if and only if the vector field 
mA is without divergence i.e. div (mA) = 0. Thus, the set of last multipliers 
is a "measure of how far away" is A from being divergence- free. 



5 



(vi) To the vector field A one may associate an adjoint vector field A*, 
acting on functions in the following manner, [21, p. 129]: 

A* (to) = —A (m) — mdivA. 

Then, another simple characterization is: m is a last multiplier for A iff m is 

a first integral of the adjoint A*. An important consequence results: the set 

of last multipliers is a subalgebra in C°° (M). 

An important structure generated by a last multiplier is given by: 
Proposition 1.4 Let m G C°° (M) be fixed. The set of vector fields 

admitting m as last multiplier is a Lie subalgebra in X (M). 

Proof Let X and Y be vector fields with the required property. Since 

[12, p. 123]: 

div [X, Y] = X (divY) - Y (divX) 

one has: 

[X, Y] (m)+mdiv [X, Y] = (X (Y (m)) + mX {divY))-{Y (X (to)) + mY (divX)) = 
= {-divY ■ X (m))-{-divX ■ Y (to)) = divY -mdivX-divX-mdivY = 0. □ 

2 Last multipliers on Poisson manifolds 

Let us assume that M is endowed with a Poisson bracket {, }. Let / G 
C°° (M) and Af e X (M) be the associated Hamiltonian vector field of the 
Hamiltonian f cf. [12]. Recall that given the volume form V there exists a 
unique vector field X v , called the modular vector field, so that ([9], [24]): 

div v A f = X v (/) . 

The triple (M, {, }, V) is called ([24]) unimodular if Xy is a Hamiltonian 
vector field, A p of p G (M). 
From (1.2) there results: 

= A f (m) + mX v (/) = -A m (/) + mAV (/) 

which means: 

Proposition 2.1 m is a last multiplier of Af if and only if f is a first 
integral for the vector field mXy — A rn . In the unimodular case, m is a last 
multiplier for Af if and only if m{p, /} = {m, /}. 
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Since / is a first integral of Af we get: 

Corollary 2.2 / is a last multiplier for Af if and only if f is a first 
integral of the vector field X v . In the unimodular case, f is a last multiplier 
for Aj if and only if {p, /} = 0. 

With the Jacobi and Leibniz formulas the following consequence of Corol- 
lary 2.2 may be established: 

Corollary 2.3 Let (M, {, }) be a unimodular Poison manifold and let F 
be the set of smooth functions f that are last multipliers of Af. Then F is a 
Poisson subalgebra in (C°° (M) , •, {, }) . 

Another important consequence of Proposition 2.1 is: 

Corollary 2.4 If m is a last multiplier of Af and A g then m is a last 
multiplier of Af g . Then, if m is a last multiplier of Af then m is a last 
multiplier of Afr for every natural number r > 1. 

Let (x 1 , . . . , x n ) be a local chart on M such that V = dx 1 A ... A dx n and 

n 

the bivector n of (M, {, }) is: n = £ ir^£ A £j. Denoting tt* = £ f£ we 

i<j j=l 

have ([5, Proposition 1, p. 4]): 

n 

dx 



i=l 

and then, Proposition 2.1 and Corollary 2.2 become: 

Proposition 2.5 (i) m G C°° (M) is a last multiplier for Af if and only 

if- 

(ii) / G C°° (M) ^ a last multiplier for Af if and only if. 

df 



n 



^ n dx 1 

i=l 



0. (2.3) 



Examples: 
2.1 

After [23, p. 31] the bivector n = h(x,y)^ A ^ defines a Poisson 

•! » W 1 t limi 

c*y ' dx 

dh df dh df 
dy dx dx dy 



structure on M 2 . So, n 1 = |^,7r 2 = and then (2.3) becomes: 
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with the obvious solution f — h. Therefore, on the Poisson manifold (IR 2 , 7r) 
above the function h is a last multiplier for exactly the Hamiltonian vector 
field A h . 

2.2 Lie-Poisson structures 

Let g be an ri-dimensional Lie algebra with a fixed basis B — {ci}i<j<rt 
and let B* = {e 1 } be the dual basis on the dual Q*. Recall the definition of 
structure constants of Q: 

[ej, 6j] = Cjj-efc. 

Then, on Q* we have the so-called Lie-Poisson structure given by ([23, p. 
31]): 

tfi{x u e u )=c%x k . (2.4) 

Hence: 

n 

^-E4 ( 2 - 5 ) 

yields: 

Proposition 2.6 For a Lie-Poisson structure on Q* provided by the struc- 
tural constants (c^) a function f e C°° (C?*) «s a Zas£ multiplier for Af if and 
only if. 

± < - < 2 - 6 ' 

Particular case: n=2 

From [ei, e x ] = [e 2 , e 2 ] = 0, [ei, e 2 ] = cj 2 ei + c^ 2 e 2 the last equation reads: 

5/ _ ! 5/ _ 

For example, if c\ 2 ■ c\ 2 ^ 0, a solution of (2.7) is: 



/ = /(*i,* 2 ) = A(^ + ^) + B (2.8) 



C 12 C 12 



with A, B real constants. 
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3 The Riemannian case 



Let us suppose that a Riemannian metric g = <, > on M is given; then, there 
exists an induced volume form V g . Let wgA 1 (M) be the g-dual of A and S 
the co-derivative operator 5 : A* (M) — ► A* -1 (M). Then: 

j divv g A = —5u 
[ A (m) = g _1 (dm, uj) 

and condition (1.2) means: 

g^ 1 (dm,u)) = mSu. 

Supposing that m > it follows that m is a last multiplier if and only 
if u belongs to the kernel of the differential operator: g~ l (dlnm, •) — 5 : 
A 1 (M) -> A = C°° (M). 

Now, assume that the vector field A admits a Helmholtz type decompo- 
sition: 

A = X + Vu (3.1) 

where X is a divergence-free vector field and u G C°° (M); for example, if 
M is compact such decompositions always exist. From divv g S7u = Au, the 
Laplacian of u, and Vw (m) =< Vu, Vm > there follows that (1.2) becomes: 

I(«) + <V«,Vm>+m-A« = (3.2) 

Example 3.1 

u is a last multiplier of A = X + Vw if and only if: 

X (u) = -u ■ Au- < Vu, Vu > . 

Suppose that M is a cylinder M = I x N with I C M and N a (n - 1)- 
manifold; then, for X = — | e (/) which is divergence- free with respect 
to V — dt A Vjv with VJv a volume form on X, the previous relation yields: 

u t = 2(u- Au+ < Vu, Vu >) . 

By the well-known formula ([18, p. 55]): 

< V/, Vg >= l - (A (fg) - f ■ Ag - g ■ Af) (3.3) 
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the previous equation becomes: 



ut = A (u 2 ) (3.4) 

which is a nonlinear parabolic equation of the type of porous medium equa- 
tion([l]). 

Example 3.2 

Returning to (3.1), suppose that X = 0. The condition (3.2) reads: 

m ■ Au+ < Vu, Vm >= (3.5) 
which is equivalent, via (3.3) to: 

A {um) + m ■ Au = u • Am. (3.6) 

Adding to (3.6) a similar relation with u replaced by m leads to the 
following conclusion: 

Proposition 3.3 Let u,m G C°° (M) such that u is a last multiplier of 
Vm and m is a last multiplier of Vm. Then u ■ m is a harmonic function 
on M. u G C°° (M) is a last multiplier of A = Vm if and only if u 2 is a 
harmonic function on M. 

If M is an orientable compact manifold then, from Proposition 3.3, it 
follows that u 2 is a constant which implies that u is a constant, too. But 
then A = V« = 0. Therefore, on an orientable compact manifold, a function 
cannot be a last multiplier of its gradient vector field. 

If (M,g) = (M n , can), there are two classes (with respect to the sign ±) 
of radial functions with harmonic square: 

i) n = 2 

u± (r) = ±^fdlnr + C 2 

ii) n — 1 and n > 3 

if±(r) = ±y/C 1 r*- n + C 2 

with Ci, C 2 real constants and r = a/ (re 1 ) 2 + . . . + (x n ) 2 . 

Example 3.3. The gradient of the distance function with respect 
to a two dimensional rotationally symmetric metric 

Let M be a two dimensional manifold with local coordinates (t, 9) en- 
dowed with a rotationally symmetric metric g = dt 2 + ip 2 (t)d9 2 cf. [18, p. 
11]. Let u G C°° (M), it (t, 6 1 ) = t, which appears as a distance function with 
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respect to the given metric. Then, Vw = and Au = equation (3.5) 
is: 

ip' (t) dm 
m- ^-7^ + — = 

ip (t) dt 

with solution: m = m(t) = p' 1 (t). 

This latter function has a geometric significance: let T — T (t) be an 
integral of m i.e. ^ = m = ^y. Then, in the new coordinates (T, 6) the 
given metric is conformally Euclidean: g = p> 2 (t) (dT 2 + d6 2 ) where t — t (T). 

Example 3.4 

Another way to treat the case X = is via equation (1.4): 

div{mVu) = (3.7) 

in which the transformation v = u^fm (recall that we search for m > 0) is 
considered. From div(mSJ -j=) = 0, it results div{^mVv — vV^/m) = i.e.: 

y^mAv = vA^/m (3.8) 

which yields: 

Proposition 3.4 Let a > 0,6 be solutions of Helmholtz [particularly 
Laplace) on the Riemannian manifold (M,g). Then, m = a 2 is a last multi- 
plier for the gradient vector field of function u = K 



4 Applications to gas dynamics 

Consider again the Riemannian manifold (M, g). Set m G C°°(M ) and recall 
according to [20, p. 62]: 

Definition 4.1 A form uj is said to be: 

(i) m-coclosed if 5(mu>) = 0, 

(ii) m-harmonic if it is closed and m-coclosed. 
An important result from the cited paper is: 

Proposition 4.2 // the 1-form co is m-harmonic, then, locally, uj = d<p, 
where <p G C°°(M) satisfies: 

5(md4>) = 0. (4.1) 
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In local coordinates (x 1 , ...,x n ) on M this equation has the form: 

-^=^-(JaWgg ij m^r) = 

and for a flat M the last equation is exactly the classical gas dynamics 
equation conform [20, p. 63]. 

But (4.1) is exactly divv g {mS7 4>) = 0, which means that m is a last mul- 
tiplier for the gradient vector field V0. So, the last result can be rephrased: 

Proposition 4.3 If the 1-form uj is m-harmonic, then, locally, u = d<p, 
where <fi G C°°(M) with V0 having m as a last multiplier. 

Due to the local character of previous result for this setting, the Propo- 
sition 1.4 can be improved: 

Proposition 4.4 Let the 1-form uj be m-harmonic and a,b G C°°(M), 
such that uj = da = db. Then m is a last multiplier for Va and Vb but first 
integral for the vector field [Va, V6]. 

Proof We have: 

[Va, V6] (m) = Va (mScu) - Vb (m5u) = m [Va (Sou) - Vb (Sou)] . 
But Va (Slo) = V6 (Slo) = g- 1 (dSco, uj). □ 
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